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ABSTRACT 

The NEA 2001 SN263 is a triple system of asteroids and it is the target of the ASTER MIS¬ 
SION - First Brazilian Deep Space Mission. The announcement of this mission has motivated 
a study aimed to characterize regions of stability of the system. Araujo et al., (2012), cha¬ 
racterized the stable regions around the components of the triple system for the planar and 
prograde cases. Through numerical integrations they found that the stable regions are in two 
tiny internal zones, one of them placed very close to Alpha and another close to Beta, and 
in the external region. For a space mission aimed to place the probe in the internal region 
of the system those results do not seem to be very interesting. Therefore, knowing that the 
retrograde orbits are expected to be more stable, here we present a complementary study. We 
now considered particles orbiting the components of the system, in the internal and external 
regions, with relative inclinations between 90° < I ^ 180°, i.e., particles with retrograde 
orbits. Our goal is to characterize the stable regions of the system for retrograde orbits, and 
then detach a preferred region to place the space probe. For a space mission, the most in¬ 
teresting regions would be those that are unstable for the prograde cases, but stable for the 
retrograde cases. Such configuration provide a stable region to place the mission probe with 
a relative retrograde orbit, and, at the same time, guarantees a region free of debris since they 
are expected to have prograde orbits. We found that in fact the internal and external stable re¬ 
gions significantly increase when compared to the prograde case. For particles with e = 0 and 
/ = 180°, we found that nearly the whole region around Alpha and Beta remain stable. We 
then identified three internal regions and one external region that are very interesting to place 
the space probe. We present the stable regions found for the retrograde case and a discussion 
on those preferred regions. We also discuss the effects of resonances of the particles with 
Beta and Gamma, and the role of the Kozai mechanism in this scenario. These results help 
us understand and characterize the stability of the triple system 2001 SN263 when retrograde 
orbits are considered, and provide important parameters to the design of the ASTER mission. 
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1 INTRODUCTION 

The population of NEAs (Near-Earth Asteroids) is composed 
by asteroids orbiting the Sun with perihelion q ^ 1.3 AU. They are 
classified according to their orbital characteristics as: Amor (a > 1 
AU and 1.017 ^ q ^ 1.3 AU), Apollo (a > 1 AU and q ^ 1.017 
AU), Atens (Q > 0.983 AU) and IEO - interior to the Earth’s orbit 
(Q < 0.983 AU), where a is the semi-major axis and Q is the 
aphelion distance of the orbit of the asteroid. 

The population of NEAs is composed mostly by asteroids co¬ 
ming from the Main-Belt asteroids. According to Morbidelli et al., 
(2002), the three most important known sources of NEAs are the 
vq secular resonance with Saturn, the 3:1 mean-motion resonance 
with Jupiter, and the Mars-crossing region. 
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Once belonging to the NEAs’ population the close encoun¬ 
ters of those asteroids with the terrestrial planets are quite frequent. 
Such characteristic makes the NEAs more accessible and so, inte¬ 
resting objects for exploration through space missions. 

Example of a successful mission to a NEA is the Hayabusa 
mission, performed by the Japan Aerospace Exploration Agency 
and that explored the NEA (25413) Itokawa in September, 2005 
(Yoshikawa, 2006). The OSIRIS-REx is a mission developed by 
NASA designed to return to the Earth with a sample of the NEA 
(101955) Bennu. The mission is planned to be launched in 2016, 
and reach the NEA in 2018 (Hergenrother et al., 2014). Besides 
those examples, we can also cite the Don Quijote program (Gal¬ 
vez, 2003), the ISHTAR program (D’arrigo, 2003), the projects SI¬ 
MONE (Wells, 2003) and MARCO POLO (Amata, 2009), all of 
them designed to send space probes to a NEA. 

The multiple systems in the NEA population are especially in¬ 
teresting for a space mission since they increase the observational 
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Table 1 . Physical and orbital data of the three components of the system 2001 SN263 


Body 

Mass^ 1 ) 

Radius 

Period^ 1 ) 

a(D 

e (!) 

fi 1 ) (*) 

Alpha 

917.47 x 10 10 kg 

1.3 km 

2.80 years 




Beta 

24.04 x 10 lo kg 

0.39 km< 2 ) 

6.23 days 

16.63 km 

0.015 

0.0° 

Gamma 

9.77 x 10 10 kg 

0.29 km< 2 ) 

0.69 days 

3.80 km 

0.016 

« 14° 


1 - Fang et al., (2011). 

2 - Araujo et al., (2012) 

(*) With respect to the equator of Alpha. 


possibilities. According to Bottke & Melosh (1996a), approxima¬ 
tely 15% of the NEAs are expected to be multiple system. Curren¬ 
tly there are known 48 multiple systems in the NEAs population 1 , 
being 46 binaries, and two triple systems: 2001 SN263 (Nolan et 
al., 2008), and 1994 CC (Brozovic et al., 2009). 

The NEA 2001 SN263 is the target of the ASTER MISSION - 
First Brazilian Deep Space Mission (Sukhanov et al., 2010). It is a 
triple system with semi-major axis 1.99 AU, eccentricity 0.48 and 
orbital inclination 6.7°. The components of the system have diame¬ 
ters of about 2.6 km (Alpha), 0.8 km (Beta) and 0.6 km (Gamma). 
With respect to the major body Alpha, Beta has a semi-major axis 
of 16.633 km - period of about 6.2 days, and Gamma has a semi¬ 
major axis of 3.804 km - period of about 0.7 days (Fang et al., 
2011). See Tab. 1 for more details. 

The announcement of the ASTER MISSION has motivated 
a study aimed to characterize regions of stability of this system. 
Araujo et al., (2012) characterized the stability regions around the 
components of the triple system for the prograde case. They per¬ 
formed numerical integrations, for a time span of two years, of a 
system composed by the Sun, the planets Earth, Mars, Jupiter, the 
three components of the triple system, and thousands of particles 
randomly distributed within the internal and external regions of the 
triple system, with relative inclinations between 0° ^ I ^ 90°. 

For the planar case, it were found stable regions placed very 
close to Alpha and Beta. When the inclinations of the particles were 
considered, the effects of the Kozai mechanism became observable 
turning the whole region around the three bodies unstable for par¬ 
ticles with relative inclination 30° < I ^ 90°, excepted for region 
1, where some particles survived. It were also characterized mean 
motion resonances with Beta and Gamma in the internal region. It 
was obtained that particles with a w 1.8 km, experience the effects 
of a 3:1 resonance with Gamma, while the particles with a ^ 7.8 
km and a ^ 7.9 km, experience the effects of a 3:1 resonance with 
Beta and of a 1:3 resonance with Gamma, respectively. For the ex¬ 
ternal region of the system (besides Gamma), it was found that the 
region is predominantly stable, and that the variation of inclination 
does not affect the stability of the particles. 

For a space mission aimed to place the probe in the internal 
region of the system those results do not seem to be very interes¬ 
ting, since the stable regions found are tiny and also because those 
are the regions more likely to have debris. Therefore, knowing that 
stable regions are greater for retrograde orbits (Hamilton & Burns 
, 1991; Hamilton & Krivov , 1997), here we present a complemen¬ 
tary study, considering now particles orbiting the components of 


1 Available at http://echo.jpl.nasa.gov/~lance/binary.neas.html. This num¬ 
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the system in the internal and external regions, with inclinations 
between 90° < I ^ 180°, i.e., considering the retrograde case. 

From the astronomy point of view, the analysis made before 
(Araujo et al., 2012) for prograde orbits are more relevant because 
it is believed that the potential dust orbiting the system came from 
the original disk that formed the bodies, so they orbit in the pro¬ 
grade direction. In this case, retrograde orbits are not very common 
for natural bodies, and they are valid only for particles that are cap¬ 
tured by the system, so in a small number. On the other side, if the 
main goal is to study the system from the perspective of sending 
a spacecraft to visit the bodies, the retrograde orbits becomes an 
interesting option. The spacecraft can easily be placed in prograde 
or retrograde orbits, with about the same expenditure of fuel. The¬ 
refore, it is important to study the stability of retrograde orbits. The 
regions of stable retrograde orbits can be useful to place a space¬ 
craft in orbits that does not need major station-keeping maneuvers 
in order to observe the three bodies of the system. 

Since the mission is planned to stay around one year in ope¬ 
ration, a study of the stability of retrograde orbits are made for two 
years, that is a time that allows the designed mission with some 
extra time, in case the equipments last longer than planned. The 
idea is to find regions of stable retrograde orbits in locations where 
the prograde orbits are not stable. Those orbits allow the spacecraft 
to be in stable orbits in regions unlikely to have particles that may 
damage the spacecraft. Those regions are identified and shown in 
detail in the present paper. We present the regions of stability found 
for the retrograde case, and compare them with the regions found 
for the prograde case. Our goal is to complement the results pre¬ 
sented by Araujo et al., (2012), given a complete scenario on the 
dynamic of particles in the neighborhood of the system, and thus, 
provide important parameters that may guide the design of the AS¬ 
TER mission. 

Only the gravitational perturbations were considered in the 
present paper and in the previous one. The solar radiation pres¬ 
sure is a relevant mechanism to be considered in a more complete 
analysis since it is known that depending on the area-to-mass ra¬ 
tio it may result in significant eccentricity increase, and so, affect 
the stability. The analysis of how relevant this mechanism is for the 
specific problem of a space mission aimed to explore this triple sys¬ 
tem is the subject of a future work. Such analysis must be guided 
by the results presented here and by Araujo et al., (2012). 

The structure of this paper is as follows. In Sec. 2, we present 
the problem of the stability around the triple system 2001 SN263 
for the inclined retrograde case. We present the initial considerati¬ 
ons (Sec. 2.1), the adopted method (Sec. 2.2), and the results ob¬ 
tained for all the regions (Sections 2.3 to 2.6). In Sec.3, we com¬ 
pare the results presented here with those presented by Araujo et 
al., (2012). Based on this comparison, we also discuss the regions 
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within the triple system that are more interesting to place a space 
probe. In Sec. 4, we present the conclusions of the work. 


2 INCLINED RETROGRADE CASE 

In order to achieve our goal, we first seek characterize the sta¬ 
ble and unstable regions in the region around the components of 
the triple system 2001 SN263, for the retrograde case. We based 
the development of this work on the initial conditions and on the 
method adopted by Araujo et al., (2012), as follows. 


Region 1 


Gamma 

Alpha 



Region 2 


0 Region 

Beta 


2.1 Initial conditions 

We considered a system composed by the Sun, the planets 
Earth, Mars and Jupiter, the three components of the system 2001 
SN263, and by thousands of particles randomly distributed around 
Alpha, Beta and Gamma (internal region), and around the whole 
system (external region). 

The internal region was divided into four regions, defined by 
the approximated Hill’s radius of each body (see Fig. 1). The re¬ 
gions 1 and 2 are the regions where Alpha is gravitationally domi¬ 
nant. The region 3 is the region which Beta dominates. The limit 
of this region was calculated for the problem involving Alpha and 
Beta, resulting in a Hill’s radius of 3.4 km. The region 4 is the 
region which Gamma theoretically dominates. The limit of this re¬ 
gion was calculated for the problem involving Alpha and Gamma, 
resulting in a Hill’s radius of 0.6 km. The region 1 is limited by the 
orbit of Gamma, and the region 2 is limited by the orbit of Beta. 
The limits of the external region are better explained in Sec. 2.6. 

We randomly distributed particles within the regions 1, 2, 3, 
and in the external region. The region 4 was not considered since 
it is too small, and would become even smaller if the perturbations 
from Beta were considered. 

All the particles started with an eccentricity going from e = 
0.0 until e = 0.5, and with a random angular distribution: 0° ^ 
/ ^ 360°, 0° ^ uo ^ 360° and 0° ^ Q < 360°, where / is 
the true anomaly, u is the argument of the pericentre and Q is the 
longitude of the ascending node. The radial distribution of those 
particles depends on the region being considered. 

Regarding the inclination of the particles, we now considered 
the inclined retrograde case, meaning that all particles were distri¬ 
buted in the internal and external regions with an inclination with 
respect to the equator of Alpha given in the range 105° ^ I ^ 
180°, taken every 15°. 

We considered the oblateness of Alpha, with J 2 = 0.013 
(Fang et al., 2011). The obliquity of Alpha was also considered 
and was determined through the pole solution, presented by Fang 
etal., (2011) 

Beta and Gamma were assumed to be spherical bodies with 
radius given in Tab. 1. A non-sphericity of Gamma would not af¬ 
fect our results since we did not consider particles around this body. 
A non-sphericity of Beta may change the results for region 3, espe¬ 
cially for those particles closer to this body. Such exception must 
be taken into account in further studies or in practical applications. 

2.2 Method 

The method consisted on numerically integrate the equations 
of motion of the gravitational problem of N-bodies. The numerical 
integrations were performed with the Gauss-Radau numerical in¬ 
tegrator (Everhart, 1985), for a time span of 2 years. This interval 


Figure 1 . Representation of the components of triple system 2001 SN263 
and their regions of influence. The blue circles represent the Hill’s radii of 
Beta and Gamma. The red dotted circles represent the collision-lines with 
Gamma and Beta, and by definition, the limits of the internal regions 1 and 
2. Reproduced from Araujo et al. (2012). 

corresponds to about 100 orbital periods of Beta, and about 1000 
orbital periods of Gamma. The time step of the numerical integra¬ 
tions was of 0.04 days, for the internal and the external regions. 

The physical data for Alpha, and the orbital elements of Beta 
and Gamma are given by Fang et al., (2011), for the epoch MJD 
54509 in the equatorial frame of J2000 (see Tab. 1). The orbital 
elements of Alpha, Earth, Mars and Jupiter were obtained through 
the JPL’s Horizons system for the same epoch, also for the equato¬ 
rial frame. The planet Saturn was not considered. Due to the short 
time span of the numerical integration its gravitational perturbation, 
including the effects of the tsq secular resonance, is negligible. 

Throughout the numerical integrations, the particles could col¬ 
lide with Alpha, Beta and Gamma, or be ejected from the system. 
The collisions were defined by the physical radius of the three com¬ 
ponents. The distance of ejection d was defined for each region. For 
regions 1 and 2 it is defined by the collision-lines with Gamma and 
Beta, being d > 3.804 km for region 1 and d > 16.633 km for 
region 2. The ejection limit for particles in region 3 is given by 
the approximated Hill’s radius of Beta, being d > 3.4 km. Fig. 1 
shows those limits. The ejection distance for the external region is 
presented in Sec. 2.6. 

The stability and instability are defined by the number of par¬ 
ticles that survive throughout the numerical integration. Initial con¬ 
ditions that result in any loss of particles, by collision or ejection, 
define a stable region. Otherwise, we define the unstable regions. 

Following, we discuss special features of the method descri¬ 
bed above, and we present the results obtained for the internal re¬ 
gion (regions 1,2 and 3) and for the external region. 

2.3 Region 1 

In region 1, we considered particles orbiting Alpha with semi¬ 
major axis 1.4 ^ a ^ 3.2 km, taken every 0.2 km, and with eccen¬ 
tricity 0.0 ^ e ^ 0.5, taken every 0.05. For each combination of a 
and e, we considered 100 particles with random values of /, u and 
Q and with inclination going from 105° to 180°, taken every 15°. 
This combination of initial conditions resulted in a total of 11, 000 
particles distributed within this region, for each inclination value. 
The orbits were numerically integrated according to the method 
presented in Sec. 2.2. 

The diagrams of Fig. 2 show the stable and unstable regions 
found in region 1 for the retrograde case. They are composed by a 
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Figure 2. Diagram of stability of region 1 (around Alpha), for a time span of 2 years, a ) I = 105°, b ) I = 120°, c ) I = 135°, d ) I = 150°, e) / = 165°, 
f ) I = 180°. The yellow boxes marked with the small black point indicate the cases of 100% of survival. The white lines indicate the limits of the region. On 
the left is the collision-line with Alpha. On the right is the collision-line with Gamma, denoting the ejection distance d. 
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grid of semi-major axis versus eccentricity, with each combination 
of a and e represented by small squares. Each of those boxes hold 
the information of 100 particles that share the same initial values 
of a and e. The colors of the diagram indicate the percentage of 
survivors, going from 90% — 100% of surviving particles (yellow) 
to less than 10% of surviving particles (dark red). The small black 
points indicate that 100% of the particles survived (stability). The 
limits of the region are represented by the white lines. They indicate 
the combination of a and e such that the pericenter of the orbits of 
the particles are inside the body Alpha, meaning collision (on the 
left), or such that the apocenter of the orbit is beyond the ejection 
distance (on the right). 

In region 1, the stable region substantially increases for incli¬ 
nation values / ^ 135°. This effect is related to the Kozai me¬ 
chanism (Kozai, 1962). For the retrograde case the critical angle 
of Kozai is / « 141°. Particles with inclination values lower or 
near this critical angle suffer the action of this mechanism, which 
is known to generate oscillations of the eccentricity and of the mu¬ 
tual inclinations, leading to the noted instability. 

The effects of the Kozai resonance is entirely noted comparing 
the results presented in the diagrams of Fig. 2 with the results pre¬ 
sented by Araujo et al., (2012) for region 1. There is a symmetry for 
the prograde and retrograde cases, where particles with inclination 
values near the interval of the Kozai cycles (39, 2° ^ ^ 141°) 

presented the predicted instability. Out of this interval we have lar¬ 
ger regions of stability. However, we see that for the retrograde 
cases, the stable conditions are much more prominent, as expected. 
The instability appears only for the particles with initial conditions 
really close to the limits of the regions. 

It can also be noted in these diagrams the occurrence of a gap 
in the stable regions for particles with a = 2.4 km. Such characte¬ 
ristic suggests the action of a resonance. This was investigated and 
is presented as follows. 


2.3.1 Resonance in region 1 - retrograde case 


We search for a commensurability of mean motion between 
particles with semi-major axis near a = 2.4 km, with Beta or 
Gamma. The mean motion was calculated through (Murray and 
Dermott, 1999): 


n — 






( 1 ) 


where G is the gravitational constant, m p is the mass of the pri¬ 
mary, R p is the radius of the primary, and J 2 is the gravitational 
harmonic of the primary related with the oblateness of this body. 

We found that particles with a = 2.398 km have a 2:1 com¬ 
mensurability of mean motion with Gamma. We then search the re¬ 
sonant angle for retrograde orbits (see Morais & Giuppone (2012) 
and Morais & Namouni (2013)). 

We found that the resonant angle ip = 2X — X' — 3 vj' presents 
an intermittent behavior (circulating and librating), as can be seen 
in Figs. 3a and 3b. In this equation, A' and w' is the mean longitude 
and the longitude of pericenter of the particle, while A is the mean 
longitude of Gamma. 

Therefore, particles with semi-major axis near a = 2.4 km 
suffer the effects of a 2:1 resonance with Gamma. As a conse¬ 
quence, the particles on this region are perturbed in such way that 
the semi-major axis remains almost constant while the eccentri¬ 
city increases (Figs. 3c and 3d), and then, the particles cross the 
collision-line of this region, giving rise to the noticed instability. 


Table 2. Mean motion resonances in region 2 


Pertuber 

Semi-major axis (km) 

Order 

Resonant angles* 

Gamma 

6,036 

1 : 2 

cp = 2X' — A — 3vj' 

Gamma 

7, 909 

1 : 3 

p> — 3X! — A — A.vj 

Beta 

7, 997 

3 : 1 

ip = 3A - A' - 4£T 

Beta 

10,477 

2 : 1 

cp = 2A — X' — 3vo' 


* A' and vj' are the mean longitude and the longitude of pericenter of the 
particle. A and vj are the mean longitude and the longitude of pericenter 
of the satellites. 


2.4 Region 2 

In region 2 the particles orbit Alpha with 4.5 ^ a ^ 13.5 km, 
taken every 0.2 km, and with 0.0 < e ^ 0.5, taken every 0.05. For 
each combination of a and e, we considered 100 particles with ran¬ 
dom values of /, oj and O and with inclination going from 105° to 
180°, taken every 15°. This combination of initial conditions resul¬ 
ted in a total of 50, 600 particles distributed within this region, for 
each inclination value. The orbits of these particles were integrated 
following the method presented in Sec. 2.2. 

The stable and unstable regions found for those conditions are 
presented in the diagrams of Fig. 4. As presented in section 2.3, we 
adopted again a grid of semi-major axis versus eccentricity and the 
color code is also the same. The limits of the region are represen¬ 
ted by the white lines. They indicate the combination of a and e 
such that the particles cross the collision line with Gamma, i.e., the 
pericenter of the orbits are smaller than 3.804 km (on the left), or 
such that the apocenter of the orbit is beyond the ejection distance 
d = 16.633 km (on the right). 

We see that for the cases with / ^ 135° the region 2 pre¬ 
sents significant stable regions. The particles with I > 135° have 
already crossed the critical angle of Kozai or are near to the limit 
I « 141°. Thus, the effects of such perturbation are no longer ap¬ 
plicable, given rise to the noted stability. From Araujo et al., (2012), 
we see that, in fact, particles with 39, 2° < I ^ 141° are under the 
Kozai mechanism action. 

In contrast to the results presented by Araujo et al., (2012), 
showing that the region 2 is predominantly unstable for the planar 
and prograde cases, here, for the retrograde case, we found that 
this region is predominantly stable. For initial conditions out of the 
Kozai mechanism action, almost the entire region is stable. In fact, 
the instability appears only for the particles with initial conditions 
really close to the limits of the region (escape from the region or 
collisions due to crossing orbits with Gamma or Beta). 

In this region we identified the action of mean motion reso¬ 
nances of particles with Beta and Gamma. Tab. 2 brings the details 
of those resonances. They were found following the procedure des¬ 
cribed in Sec. 2.3.1. For all of them the resonant angles librate, as 
exemplified in Fig. 5 for the resonance 3:1 with Beta. The effects 
of those resonances on the particles placed in region 2 are noted 
through small gaps within the stable cases, which is more evident 
for higher values of eccentricities. 


2.5 Region 3 

In region 3 the particles orbit Beta with 0.8 ^ a < 3.4 km, 
taken every 0.2 km, and with 0.0 < e < 0.5, taken every 0.05. 
As previously done for the other regions, we considered 100 parti¬ 
cles with random values of /, uj and Q and with inclination going 
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Figure 3. Particle in region 1 with a = 2.4 km, e = 0, 0 e / = 150°. a) Evolution of the resonant angle tp = 2A — A' — Suj' for t = 2 years, b) Zoom 
of Figure (a) showing some of the libration regions. Time going from 1.4 to 1.55 years (« 79 orbital periods of Gamma), c) Evolution of semi-major axis, d) 
Evolution of eccentricity. 


from 105° to 180°, taken every 15°, for each combination of a 
and e. This combination of initial conditions resulted in a total of 
15,400 particles distributed within this region, for each inclination 
value. The orbits of these particles were also integrated following 
the method presented in Sec. 2.2. 

The stable and unstable regions found for those conditions are 
presented in the diagrams of Fig. 6. We adopted the same grid 
of semi-major axis versus eccentricity and color code for the di¬ 
agrams. The limits of the region are represented by the white lines. 
They indicate the combination of a and e such that the pericenter 
of the orbits of the particles are inside the body Beta meaning col¬ 
lision (on the left), or such that the apocenter of the orbit is beyond 
the ejection distance d = 3.4 km for region 3 (on the right). 

Similar to the results for regions 1 and 2, we note that the 
stable regions appear only when the relative inclination is about to 
reach the critical angle ofKozai I « 141°. 

The effect of the Kozai mechanism is more prominent for the 
particles in this region. From the diagram of Fig. 6 and from the 
diagrams presented by Araujo et al., (2012), for the region 3 and 
for the prograde case, we see that the stable regions completely 
vanish for inclination values near the interval 39, 2° ^ I ^ 141°. 

Similar to the results obtained for regions 1 an 2, we also note 
a substantial increase of the stable region in region 3, when com¬ 


paring the retrograde case with the prograde case (Araujo et al., 
2012). As found for regions 1 and 2, here we also found that the 
instability appears only for initial conditions really close to the li¬ 
mits of the region. Thus, for inclination values out of the Kozai 
mechanism action, we found that retrograde orbits in region 3 are 
predominantly stable. 


2.6 External Region 

We extend the analysis of the stability of the triple system 
2001 SN263 to the external region. 

The lower limit of the external region is d = 20 km from 
Alpha. This distance corresponds approximately to the distance 
Alpha-Beta added to the Hill’s radius of Beta. The upper limit was 
determined taking into account the approximated Hill’s radius of 
the whole system. Considering a body with mass equal to the sum 
of the mass of Alpha, Beta and Gamma, and the heliocentric orbit 
of the system, Araujo et al., (2012) calculated that the Hill’s radius 
of such body is Rhui ~ 180.0 km at the perihelion of the orbit 
and Rhui ~ 500.0 km at the aphelion. 

According to Domingos et al., (2006) the predicted limit of 
stability for retrograde orbits is of about one Hill’s radius. Thus, 
the upper limit of the external region was chosen to be d = 180.0 
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Figure 4. Diagram of stability of region 2 (between Alpha and Beta), for a time span of 2 years for / = 105°, / = 120°, / = 135°, / = 150°, / = 165° 
and / = 180°, from top to bottom. The white lines indicate the limits of the region.The yellow boxes marked with the small black point indicate the cases of 
100% of survival. 


km (1 Rhui at perihelion, when the system is more perturbed). 
Beyond this distance the particles are considered ejected from the 
system. 

Thus, in the external region we have particles orbiting Alpha 
with 20.0 ^ a ^ 180.0 km, taken every 1.0 km, and with 
0.0 ^ e ^ 0.5, taken every 0.05. As previously done, for each 


combination of a and e, we considered 100 particles with random 
values of /, uo and Q and with inclination going from 105° to 180°, 
taken every 15°. This combination of initial conditions resulted in 
a total of 177,100 particles distributed within this region, for each 
inclination value. The orbits of these particles were integrated fol¬ 
lowing the method presented in Sec. 2.2. 
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Figure 5. Particle in region 2 with a = 7.997 km, e = 0, 0 e I = 135°. a) Evolution of the resonant angle tp = 3A — A' — 4 Q! for t = 2 years, b) Zoom of 
Figure (a) showing some of the libration regions. Time going from 0.72 to 1.0 years (« 17 orbital periods of Beta). 


The stable and unstable regions found for the external region 
are presented in the diagrams of Fig. 7. The limits of the region 
are represented by the white lines. They indicate the combination 
of a and e such that the particles cross the collision line with Beta, 
i.e., the pericenter of the orbits are smaller than 16.633 km (on the 
left), or such that the apocenter of the orbit is beyond the ejection 
distance d = 180 km (on the right). Comparing those diagrams 
with each other, we note that there is no significant difference in the 
stable and unstable regions. The instability is observable only in the 
border of the region. As previously found in Araujo et al., (2012) 
for the planar and prograde cases, we concluded that the external 
region of the system is predominantly stable for the retrograde case, 
and that the stability of the external region of the triple system is 
not affected by the inclination variation. 


3 PROGRADE AND RETROGRADE CASES IN A SPACE 
MISSION SCENARIO 

For a space mission, the most interesting regions would be 
those that are unstable for the prograde cases, but stable for the re¬ 
trograde cases. Such configuration provide a stable region to place 
the mission probe for a period of 2 years with a relative retrograde 
orbit. At the same time, if we take into account that possible de¬ 
bris located in the neighborhood of the components of the system 
are more likely to have prograde orbits, we may have a region free 
of debris since they are expected to be quickly removed from the 
system due to the instability. 

The diagrams in Fig. 8 help us to compare the results for pro¬ 
grade and retrograde cases. They were done considering only the 
particles in the internal regions of the triple system with e = 0.0, 
and with 7 = 0° (Fig.8a), and with I = 180° (Fig.8b). These dia¬ 
grams represent the radial distribution of the particles belonging to 
the set of initial conditions (a x e) that after the numerical integra¬ 
tion resulted in 100% of survival (yellow region marked with the 
small black dots). 

For the planar prograde case we may consider the Hill stabi¬ 
lity, i.e., a stability criterion that consider the restricted three-body 
problem and that provides a formula to calculate A - the mutual 
orbital separation - such that the orbits remain stable (Gladman, 
1993). Being Alpha the central body, we calculate A considering 


the problems: Alpha-Beta-Particle (/ii = mp/rria) and Alpha- 
Gamma-Particle (fjbi — m 7 /m a ), where the particle have the ex¬ 
ternal orbit and A « 2.4(/xi) 1//3 . For the Alpha-Beta-Particle pro¬ 
blem, we found that A « 11.8 km, meaning that particles beyond 
Beta must have stable orbits for mutual distance Beta-particle gre¬ 
ater than this value. For the Alpha-Gamma-Particle problem, we 
found that A ^ 2.0 km. Although this analytical approach does 
not take into account a third disturbing body and the effects of res- 
sonances, we found that the results for the prograde case are in 
agreement with these predictions, as can be seen in the diagram of 
Fig.8a. We see a tiny stable zone located beyond Gamma starting 
at 9.1km (A « 5.3 km), and beyond Beta the stable orbits appear 
for particles with orbital radius of 35 km (A « 18.4 km). 

From the diagram in Fig.8b it became clear that for the retro¬ 
grade case the whole region around the triple system 2001 SN263 
is essentially stable. Differently from the results obtained for the 
prograde case, it is clear here, for the retrograde cases, that the ins¬ 
tability appears only for those particles orbiting Alpha close to the 
orbits of Gamma and Beta. Although it is also possible to have sta¬ 
ble particles around Beta itself. 

In the diagram of Fig. 8c, we detached the preferred region 
with the orange color. We note that in regions 1 and 3 there are 
tiny zones around Alpha and Beta that meets the requirement of 
being stable for retrograde orbits and unstable for the prograde 
ones. Around Alpha, we have a region with radial length of 0,6 
km (region 1), and of 1,0 km around Beta (region 3). The region 2 
appears as the most interesting region. It was characterized as uns¬ 
table for prograde orbits, but is stable when the retrograde orbits 
are considered. Thus, we have a region with radial length of 8.2 km 
that is stable for the space probe for a period of 2 years and that 
must be free of debris. 

Therefore, we identify three internal regions that are very inte¬ 
resting to place the space probe, since they are stable for retrograde 
orbits but unstable for prograde orbits and allow the spacecraft to 
be in stable orbits close to the components of the system. The first 
one is around the main body, in region 1. It is very useful to place 
the spacecraft to observe the main body more closely. The second 
one is around the second larger body (Beta). It is good to locate the 
spacecraft if the goal is to better observe this body. The third one 
is in the regions between the two satellite bodies (region 2). It is 
a very interesting region to keep the spacecraft to observe the th- 
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Figure 6. Diagram of stability of region 3 (around Beta), for a time span of 2 years, a) / = 105°, b ) I = 120°, c ) I = 135°, d ) I = 150°, e) / = 165°, 
f) / = 180°. The white and blue lines indicate the limits of the region.The yellow boxes marked with the small black point indicate the cases of 100% of 
survival. 
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Figure 7. Diagram of stability of the external region (beyond Beta) for a time span of 2 years for I = 105°, I = 120°, I = 135°, I = 150°, I = 165°, and 
I = 180°, from top to bottom. The yellow boxes marked with the small black point indicate the cases of 100% of survival. 


ree bodies of the system. An elliptical orbit with periapsis near the 
lower limit of this region and apoapsis near the upper limit of this 
region would allow the spacecraft to pass closer to both secondary 
bodies in a single stable orbit that stays inside a region that is free 
of concentration of particles. 

The external region are essentially stable for both, prograde 
and retrograde orbits. Nevertheless, we can detach a region of about 
7 km, going from a — 27 km to a — 34 km, beyond Beta, that meet 
the preferred region requirement. This would be a very interesting 
region to place a space probe. There we may have a clean area with 
stable retrograde orbits around Alpha that allow a close observation 
of the three components of the system. 


4 CONCLUSIONS 

In the present paper we analyze stability regions of the triple 
system 2001 SN263, for retrograde orbits in the internal and exter¬ 
nal regions of the system. Based on those results, we detached the 


preferred regions, i.e, internal and external regions that are interes¬ 
ting to place a space probe aimed to explore this system. 

We applied the method adopted by Araujo et al., (2012), but 
now considering particles with retrograde orbits, i.e., particles with 
90° < I < 180° relative to the equator of Alpha, within and 
around the entire system. We analyzed the behavior of such par¬ 
ticles to determine the stable regions for the retrograde case, and 
we compared our results with the planar and prograde cases 

Comparing our results with the planar and prograde cases we 
noted, as expected (Hamilton & Bums , 1991; Hamilton & Krivov 
, 1997), a substantial increase of the stable region for both, internal 
and external regions when the retrograde orbits were considered. 

We found that nearly the whole internal region is stable, espe¬ 
cially for low eccentricities values. The instability appears only in 
the border of the regions, or when the inclination of the particles 
are such that they are under the Kozai mechanism action. The most 
prominent result was found for region 2. Comparing the simulati¬ 
ons for particles with 7 = 0° with those with I = 180°, we see that 
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(C) 


Figure 8. Comparison of the stable (yellow) and unstable (red) regions at 
the internal and external regions of the triple system 2001 SN263 for the 
circular case, a) Prograde case, reproduced from Araujo et al. (2012). b) 
Retrograde case, c) The difference between prograde and retrograde cases - 
preferred regions. 

the region that is essentially unstable for the prograde case, became 
practically fully stable for the retrograde case. 

We found that the external region is also predominantly sta¬ 
ble. The instability appears only in the borders of the region. The 
variation of inclination does not affect the stability of the particles. 

Our results showed that the retrograde orbits are an interesting 
option for the ASTER mission. We detached at least four regions 
(three internal and one external) that can be very useful to this mis¬ 
sion. All of them allow the space probe to closely observe the com¬ 
ponents of the system in stable retrograde orbits situated in regions 
expected to be free of debris. Those results provide important in¬ 
formation about the stability around the components of the system 
and substantially contribute to the design of the ASTER mission. 
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